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We present a real-time diagrammatic formalism to study adiabatic pumping in chains of tunnel- 
coupled metallic islands. This approach is based on an expansion to linear order in the frequency of 
the time-dependent parameters and on a systematic perturbation expansion in the tunnel-coupling 
strength. We apply our formalism to single-island and double-island systems. In the single-island 
setup, we find that the first-order contribution in the tunnel-coupling strength is purely due to 
the renormalization of the charging-energy gap. In the double-island system, we investigate the 
transition between weak and strong pumping. 
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I. INTRODUCTION 



The most common method to drive a current through 
a conductor is by applying a voltage. In nanostructures, 
due to the ease with which parameters can be tuned, 
another possibility opens up: pumping. In an electron 
pump a dc current can be generated in a nanoscale con- 
ductor, even at zero-bias voltage, by periodically chang- 
ing some of its properties in time. If the time depen- 
dence of the parameters is slow compared to the inter- 
nal time scales of the system, such as the dwell time 
of carriers, pumping is adiabatic. In this case, the 
pumped charge does not depend on the detailed time 
evolution of the pumping cycle but only on its area. 
Most of the theoretical efforts on pumping have been 
addressed to systems where the electron-electron interac- 
tion is weak and can, therefore, be neglected or treated 
in a mean-field approach^^^ For these weakly in- 
teracting nanoscale conductors, a well-established for- 
malism for pumping exists based on the dynamical- 
scattering approach to mesoscopic transport i£ Re- 
cently, several works have investigated pumping through 
interacting systems Most of 

these works have focused on few-electron quantum 
dot a 10 i 12 ' 13 i 15 ' 16 i 17 ' 18 since these systems are the proto- 
type for nanostructures with strong interaction effects. 
An open question regards how the pumping characteris- 
tics of an interacting quantum dot depends on the dot's 
single-particle spectrum. In order to address this prob- 
lem, we have decided to study a quantum dot with a 
dense spectrum, i.e., with vanishing level spacing. Due 
to its finite density of states, such a quantum dot is re- 
ferred to as a metallic island. Indeed, pumping in systems 
of metallic islands has been investigated experimentally 
already in the early 1990si 19 ' 20 Pothicr et al£2- realized a 
single-electron pump in a setup with two tunnel- coupled 
metallic islands in series. Pumping was obtained by 
changing in time the gate voltages of the two islands. The 
pumping cycle was chosen to enclose a degeneracy point, 
where three charge states of the double-island system are 
degenerate. This resulted in pumping of exactly one elec- 



tron per cycle. To describe this experiment, it was as- 
sumed that the system relaxes always to the ground state 
and that this transition occurs via first-order tunneling. 
Although this theoretical explanation is perfectly suited 
to describe the pumping cycle of Ref. [20|, it would fail 
in describing weak pumping or pumping cycles nonencir- 
cling a degeneracy point. 

More experiments on charge pumping through nanos- 
tructures were performed on metallic multijunction sys- 
tems^ double quantum dots in InAs quantum wires;22 
or making use of surface acoustic waves i 23 i 24 ' 25 Also spin 
pumping through a quantum dot was reported^ 

In the present paper, we develop a diagrammatic real- 
time approach to adiabatic pumping in systems consist- 
ing of metallic islands. To this aim, we follow the lines of 
Ref. [l5|, where adiabatic pumping through a single-level 
quantum dot was considered. Our approach relics on a 
systematic pcrturbativc expansion in the tunnel- coupling 
strengths and therefore is restricted to weak tunnel cou- 
pling. In particular, it is not suitable to study the very- 
low-temperature regime, where high-order tunneling pro- 
cesses become relevant. However, it treats the Coulomb 
interaction on the island exactly, allowing us to investi- 
gate effects due to strong interaction. 

First, we apply our formalism to a single- island sys- 
tem consisting of one metallic island with Coulomb in- 
teraction, tunnel coupled to two noninteracting leads. 
We compute the pumped charge up to first order in the 
tunnel-coupling strength for any pair of pumping param- 
eters chosen among the charging-energy gap of the sys- 
tem and the left and the right tunnel- coupling strengths. 
As a result, we find that the contribution in first or- 
der is due solely to the renormalization of the charging- 
energy gap. For the case of pumping with the two tunnel- 
coupling strengths, this term becomes the dominant one. 
Furthermore, we consider pumping with the charging- 
energy gaps in a more complex system consisting of two 
tunnel-coupled metallic islands called the double-island 
system. We calculate the pumped charge through this 
system in lowest order in the tunnel- coupling strength 
and investigate the issue of pumped-charge quantization. 
With our technique we arc able to describe the transition 
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from weak to strong pumping, finding charge quantiza- 
tion for strong pumping. 

This paper is structured as follows. In Sec. II 
we present the model. The technique to compute the 
pumped current is introduced in Sec. III. The results for 
single-island and double-island systems are presented in 
Sees. IV and V, respectively. The conclusions are given 
in Sec. VI. Details on the calculation of the diagrams can 
be found in Appendixes A and B. 

To keep all formulas transparent we set h = 1 through- 
out the paper. 



II. MODEL 

We consider a chain of M tunnel-coupled metallic is- 
lands. The leftmost and rightmost islands are tunnel 
coupled to metallic leads. The different metallic re- 
gions are labeled by the index r, where r = and 
r = M+l correspond to the noninteracting metallic leads 
and r = 1, . . . , M to the metallic islands. The metallic is- 
lands, for typical experimental realizations, have a dense 
single-particle spectrum and can be treated as equilib- 
rium reservoirs of fermionic degrees of freedom, described 
by a Fermi function. Each island is capacitively coupled 
to a gate voltage, which can be used as a pumping param- 
eter. The state of the chain is described by the number 
of excess charges N r on each island, |x) = |iVx, ■ ■ ■ , Nm}- 
The Hamiltonian of the system is H = Hq + H tunn . The 
operator Hq describes the system in the absence of tunnel 
coupling 



island r, i.e., [ip r ,N r ] 



Therefore. e ±?t ^ 



M+l 



<~ki 



H, 



ch • 



(l) 



=0 fc, 



where e^ nr are the single-particle energies, labeled by k 
and by the transverse-channel index n — 1, . . . , A c h an n- 
The spin degree of freedom is included in the sum over 
transverse channels. In the following, we will consider 
the many-channel case, A c haim 3> 1. The Coulomb inter- 
action is taken into account by H c \ ll which describes the 
electrostatic energy of a given charge state \x}~ This 
charging energy depends on the gate voltages applied to 
each single island. The tunnel coupling is described by 
the tunneling Hamiltonian 



M 

-fftunn = 

r=0 k,p,n 



V r (t)c{ nr c pnr+1 e -i(w +1 -*v) + H.c.1 (2) 



with tpo = (p M+i = 0. For simplicity in Eq. ^ we have 
assumed that the tunnel-matrix elements V r for tunnel- 
ing through barrier r do not depend on the quantum 
numbers k and p or on the transverse-channel index n. 
Notice that the tunneling Hamiltonian Eq. @ conserves 
the transverse-channel index. In Eq. ([2]) we have ex- 
plicitly indicated that the tunnel-matrix elements can be 
time dependent. The operator ip r is the canonical conju- 
gate of the number operator N r for the excess charges on 



are simply 

ladder operators for the charge on island r. Finally, we 
characterize the junction r between the two consecutive 
islands r and r + 1 by the tunnel-coupling strength 



«5(t) 



4n 2 Rr(t) 



(3) 



where = h/e 2 is the quantum resistance, R r the tun- 
nel resistance of the junction, A™(0) and A" +1 (0) the 
densities of states at the Fermi energy for channel n 
of leads or islands r and r + 1, respectively. We de- 
fine also the two-time coupling strength as ctQ(t,t') = 
£ n AT r "(0) N? +1 (0)V r (t)V r *(t'). 



III. FORMALISM 

We are interested in the dynamics of the excess charge 
on each island and not in the large number of nonin- 
teracting fermionic degrees of freedom which act effec- 
tively as bathsi 28 i 29 ' 30 Therefore, we trace over the states 
of the fermionic sector and we obtain the reduced den- 
sity matrix g Te d, which just describes the charge de- 
grees of freedom. In the basis defined by the states 
|x) = \Ni, . . . , Nm) (with excess charge N r on the is- 
land r), the dynamics of the diagonal and off-diagonal 
matrix elements of the reduced density matrix are un- 
coupled. Furthermore, the off-diagonal matrix elements 
do not contribute to the current. Hence, we need only to 
consider the diagonal matrix elements, which are the oc- 
cupation probabilities for the respective states and which 
we denote by p x (t) = (xlft-ed(i)lx)- The probabilities 
obey the generalized master equation 



t 

j t p(t)= J dt'W(t,t')p(t') , 



(4) 



where p is the vector of the probabilities p x and the ma- 
trix element W xx > (t, t') of the kernel W (t, t') includes all 
transitions from state x' at time t' to state x at time t. 



A. Adiabatic expansion 

We now proceed by performing an adiabatic expansion, 
i.e., a perturbation expansion in the frequency f2 of the 
time-dependent pumping fields Xi(r). To describe the 
adiabatic pumping regime it is enough to retain terms 
linear in f2. The first step consists in performing a Tay- 
lor expansion of p(t') on the right-hand side (r.h.s.) of 
Eq. ^ around t up to linear order, 



dt'W (t, t') 



(f-t) 5 P(«) 



(5) 
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This expansion is justified if the period of the pumping 
fields is much longer than the characteristic memory time 
of the system. Second, we make an adiabatic expansion of 
the kernels appearing in the generalized master equation 
as follows: 



W(t, t') -> W t (i) (t - t') + W ( t a> (t - t') 



00/ 



(6) 



The instantaneous term, (t - 
by replacing Xi(r) with Xi(r 



t'), is simply obtained 
= t). On the other 

hand, the first adiabatic correction, Wj (t — t') is ob- 
tained by performing the Taylor expansion Xi(r) — > 
Xi{t) + (r — t) £pXi(r)\ and retaining terms up to 
first order in the time derivatives. The index t in Eq. (|6|) 
indicates the time with respect to which the adiabatic 
expansion has been performed. The adiabatic expansion 
for the kernels is valid if the period of the pumping fields 
is much longer than the characteristic response times of 
the system, i.e., if <C ao max{A r , /cb? 1 }, where A r are 
the charging-energy gaps of the islands. 
Finally, we expand also the probabilities, 



complemented by the normalization conditions 
e TpO> ) _ i anc j eTpO' 1 ) _ w here e is the vec- 
tor with all components equal to one. The first adiabatic 
correction to the probabilities obey the equations 



dt 
d_ 
Tt 



(i.O) 
Pi ' 



0,1) 
Pi 



txrO, 1 ) (a,-l) 



(10) 



w 



(i,l)„(a,0) 



w 



(i,2) (a,-l) 



complemented 

l 



T ( a i 

e Pt 



W (a4) p (i,0) +aw 



by the 



4 dt Pt ' [ ' 



and e T pj a ' C 



normalization conditions 
= 0. The orders in qq 
on the r.h.s. of Eqs. (p~0|) and (fTT|) have been chosen to 
match those on the left-hand side. We notice that the 
expansion in ao starts with the order — 1. Although this 
result might seem surprising at first, it is in line with our 
expansion since p| a ' cx Q/(cto max {A r , k^T}) <C 1. 



P(*) 



.(0 



fa) 



(7) 



C. Pumped current and charge 



Notice that by construction are the stationary prob- 
abilities for the time-independent problem where all pa- 
rameters are frozen at their values at time t. Since all 
expressions now depend on the time difference t — t' it 
is useful to introduce the Laplace transform F (z) = 
J**^ dt' e" 2 ^-*') F(t - t'). In the following we 



(i/a) 



make use of the notation W._ 
J*^ dt' Wf /a) (t - t') and dW[ i] 



= W« i/a) (z) 
& W« (z) 



z = 0a 



z=0 4 



B. Perturbative expansion in the tunnel coupling 

Even after having performed the adiabatic expansion, 
the problem cannot be solved exactly due to the presence 
of the Coulomb interaction and of the tunnel couplings. 
Therefore, in addition to the expansion in the pumping 
frequency f2, we perform a perturbative expansion in the 
tunnel-coupling strength a®. This limits the range of va- 
lidity of the method to the weak tunnel- coupling regime. 

The order in the expansion in will be indicated by 
an integer superscript. For example, wj 3 ' 1 ' indicates the 
first order in ao of the adiabatic correction of the kernel. 

The perturbative expansion for the instantaneous 
probabilities is straightforward since it corresponds to the 
time-independent problem with all parameters frozen at 
time t. The instantaneous probabilities pj 1 ' ^ and p\ 
are determined from 



- wC^Ui-O) 



W 

w 



(i,2)ti>0) 
Pt 



w 



(i,l) 

Pt 



(8) 
(9) 



The pumped current through the left lead can be writ- 
ten as 



7 L (t) = e j dt'e T W L (t, t') p (t 1 ) 

J — OO 



(12) 



with e being the electron charge and W L (£,t') = 
EpP wLp (*, *') ■ Th e kernel W rp (t, t') describes all pro- 
cesses where the number of electrons entering reservoir r 
minus the ones leaving it equals p. 

We proceed for the current in the same way as for the 
master equation (0J performing an adiabatic expansion 
in the frequency O and a perturbative expansion in the 
tunnel coupling ao. Since there is no applied transport 
voltage, the instantaneous current vanishes. The adia- 
batic corrections to the current in zeroth and first order 
in ao read 



,(a,0) 



(t) 



e e T W, L(i ' 

r L(i,l) (a,0) 



pi 
V 



(13) 



W 



L(i,2) fa -1) 



Pt 



W L(M) p (U» 



aW L(i,l)rfPt 

* dt 



0\o) 



.(14) 



The lowest-order contribution to the adiabatic current, 
given by Eq. (|T3"|l . scales with fl and is independent of 
ao, in contrast to the dc current through a system with 
an applied transport voltage, which scales with ao- Since 
Q, <C «o max {A r , fceT}, the pumped current goes to zero 
for vanishing tunnel coupling as it should. The pumped 
charge Q per cycle is simply given by the integral of the 
pumped current over one period T = 



4 



-iz 





FIG. 1: Example of diagrams for the single-island system in 
first order in the tunnel coupling which are needed to calcu- 
late the kernels wf /a,1> (2) and dWf' 1} (z). The additional 
frequency lines needed for the adiabatic corrections to the 
kernels are shown as dotted lines. 



D. Diagrammatic rules 

1. Rules to calculate the instantaneous kernels 

We start with summarizing the diagrammatic rules for 
the kernel wf'™' (z) , where n indicates the order of the 
expansion in ao- Examples of first-order diagrams for the 
single-island system are drawn in Fig. [1] 

(1) Draw all topologically different diagrams with n di- 
rected tunneling lines connecting pairs of vertices 
containing lead-electron operators. Assign a reser- 
voir index r and energy cu to each of these lines. 
Assign charge states \ an d the corresponding ener- 
gies E x (t) to each element of the Keldysh contour 
connecting two vertices. Furthermore, draw an ex- 
ternal line from the upper leftmost beginning of an 
island propagator to the upper rightmost end of 
an island propagator that carries the (imaginary) 
energy — iz. 

(2) For each time segment between two adjacent ver- 
tices (independent of whether they are on the same 
or on opposite branches of the Keldysh contour) 



write a resolvent 



AE(t) 



where AE(t) is the dif- 



ference of left-going minus right-going energies (in- 
cluding energies of tunneling lines and the external 
line - the positive imaginary part of iz will keep all 
resolvents regularized) . 

(3) Each vertex containing an island operator cj^ r 

gives rise to a matrix element (x'l c in r lx)> where 
X {x') is the island state entering (leaving) the ver- 
tex with respect to the Keldysh contour. 

(4) The contribution of a tunneling line of reservoir r 
is a+ (t, u>) if the line is going backward with re- 
spect to the closed time path and a~(t,to) if it is 
going forward with <xfr(t,u)) = ±ao(t)- 
and (3 = l/{k B T). 



'cxp(±/3^ 



(5) The overall prefactor is given by (— l) b+c , 
where b is the total number of vertices on the back- 



ward propagator and c the number of crossings of 
tunneling lines. 

(6) Integrate over all the energies of the tunneling lines 
and sum over the reservoir indices. 



2. Rules to calculate the adiabatic kernels 

The possible time-dependent parameters are the 
tunnel- matrix elements V r (ti) and the charging-energy 
gaps A r (t). The matrix elements V r (ti) enter in the con- 
tractions of tunnel vertices as a r (ti, tj) oc V r (ti) V* (tj). 
The charging-energy gaps A r (t) enter in the isolated- 
island propagator exp(— i f. 3 dt'E x (<')). 

While for the instantaneous kernels all parameters were 
taken at time t, we now perform an expansion around the 
same time t and keep all contributions linear in a time 
derivative of the pumping parameters, 



al(ti,tj) ps a r (t) 



(15) 



-ifdt' E x (t') 



ti — t d«n / s ti — t dajl , . 



-i E x {t)-{t 3 -U) 



(16) 



1-i 



. fa - tf - (U - tf dE x 



dt 



(t) 



The linear factors (ti — t) appearing in these expansions 
can be included in the diagrammatic rules by introducing 
an additional external frequency line with the imaginary 
energy —izi from the vertex at t% to the rightmost vertex 
at t, performing the first derivative with respect to z,-, 
and then setting z% = 0+. Similarly, the quadratic terms 
(ti — t) 2 can be included in the diagrammatic rules by in- 
troducing an additional external frequency line with the 
imaginary energy —izi from the beginning of an island- 
propagator line at ti to the rightmost upper end of an 
island-propagator line at t, performing the second deriva- 
tive with respect to Zi and then setting zi = 0+ . The ex- 
ternal frequency lines are drawn as dotted lines in Fig. [T] 
The rules to compute the contribution to the adiabatic 
corrections W^'™" 1 due to the time dependence of cvjj(t) 
read: 

(7a) Add to all diagrams needed for (z) additional 

external frequency lines between any vertex ti and 
the upper right corner of the diagram and assign 
to them an (imaginary) energy —izi. Note that an 
external frequency line between two right corners 
of a diagram does not contribute and can always 
be omitted. 



(7b) Follow the rules (1) to (6) taking into account the 
extra lines. 

(7c) Perform a first derivative with respect to Zj and 

Sum all the contribu- 



multiply it by 



tions obtained in this way. 
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(7d) Set all the external frequencies zi and z to 0_|_. 

The contribution to the adiabatic correction W^ a ' n ^ 
due to the time dependence of the charging energy can 
be computed in a similar way: 

(8a) In addition to the external frequency lines added 
according to rule (7a), put one more external fre- 
quency line from the left corner of the diagram with 
no vertex to the right corner. 



(8b) Follow the rules (1) to (6) taking into account the 
extra lines. 



ic) Perform a second derivative with respect to Zi and 

multiply it by -| ^ (t), where \ (xO is the 
island state entering (leaving) the vertex of the 
external frequency line at with respect to the 



Keldysh contour. The term (^f^- 



is omit- 



dt \ dt 

ted if the segment associated with E x {E x i) does 
not belong to the diagram. Sum all the contribu- 
tions obtained in this way. 

(8d) Set the external frequencies Zi and z to 0+. 

Finally to compute the current, we need the matrix el- 
ements of W[ = J2 P P W[ p . For each contribution W[ p 
all diagrams for which the number of tunneling lines with 
reservoir index r running from the upper to the lower 
propagator minus the number of those with reservoir in- 
dex r running in the opposite direction equals p need to 
be added. 



IV. ADIABATIC PUMPING THROUGH A 
SINGLE-ISLAND SYSTEM 



In this section, we present the results for a single-island 
system consisting of one metallic island with Coulomb 
interaction, tunnel coupled to two nonintcracting leads 
kept at the same chemical potential. The Coulomb inter- 
action on the island is described by the charging energy 
E ch (N,N x (t)) = E C {N - N x (t)) 2 , where N x (t) denotes 
the charges externally induced via the gate voltage Vq 
and TV indicates the charge on the island. The energy 

2 

scale for the charging energy is Eq = ^ with C being 
the total island capacitance. 

At low temperatures only two different charge states 
contribute to the transport. We indicate these states 
as |0) and |1) and their energies as Eq and Ei, respec- 
tively. The occupation probabilities for these states are 
labeled po and p\ . The charging-energy gap of the island 
is defined as the energy difference A(t) = Ei{t) — Eo(t) . 
An energy sketch of the single-island system is shown in 
Fig.H 





U(t) 




c&t) 




o&t) 



lead L 



metallic 
island 



leadR 



FIG. 2: Sketch of the single-island system. The possible 
pumping parameters are the charging-energy gap A of the 
island and the tunnel-coupling strengths Oq. 



A. Pumped current in zeroth order in ao 

To calculate the pumped current in zeroth order in 
ao we need the instantaneous occupation probabilities 
in zeroth order in ao and their adiabatic correction in 
minus first order. The instantaneous probabilities can be 
obtained from Eq. Q and are simply given by Boltzmann 
factors, 



(17) 



where (3 = is the inverse temperature, E x the energy 

of the island state x, and Z = J2 X e~@ Ex the partition 
function. The instantaneous average occupation number 
of the island in lowest order in ao can be calculated from 
the instantaneous probabilities and reads 



(n) m = / + (A) 



(18) 



The adiabatic correction to the probabilities can be 
obtained from Eqs. (fTU|) and (jTTJ) and is given by 



P/ 



(a, 



2?ra(A) dt 



dim 



(19) 



where we used the notation a(cu) — a + (uj) + a~ (u>) and 
a ± (u;) = X)r a r t ( to ')- Furthermore, we define ao = a^ + 
aig". Since this probability vector is proportional to the 

time derivative of the instantaneous probabilities p[ 1,0 \ 

it is an eigenvector of the kernel Wj 1 . To compute 
the probabilities Eqs. (fl7|) and (fT9|) . we need to evaluate 
only the instantaneous kernels in first order in ao. The 
expressions for the matrix elements of are given 

in Appendix [5] . 

By inserting the results for the probabilities in 
Eq. (|13|) . it yields for the adiabatically pumped current 
through the single- island system in zeroth order ao, 



r (a,0) 



ao at 



v(i,0) 



(20) 



The average occupation number of the island in zeroth 
order ao depends only on the charging-energy gap A 
and not on the tunnel-coupling strengths. Therefore, the 
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pumped current in zeroth order ao is nonvanishing only 
if the charging-energy gap is one of the pumping parame- 
ters. The result of Eq. (|20|) has a simple interpretation: if 
the charging-energy gap is varied in time, the occupation 
of the island adjusts to the new situation and therefore 
electrons need to enter or leave the island. The fraction 
of the resulting electrical current which flows through 
barrier r is given by the ratio oi^/ao. Therefore, this 
pumping mechanism can be considered as peristaltic. 



B. Pumped current in first order in ao 



In order to compute the first-order correction to the 
pumped current, we need, besides p^ 1 ' 1 ^, also pj a '°\ 
which can be computed by means of Eq. (fTTj) and reads 



(a,0) 

Pi 



(a,0) 
Pi 



>(A) 



a(A) 



2^K[g(A)] d A a(A) - [/+(A)] 



5S"WA)]^(A) 



(25) 



We start by giving the expression for the instantaneous 
average occupation number of the island in first order in 
ao as follows: 



(n) 



0,1) _ 



a { Iq(A)} — (n 



0,o) 



v0,0) 



- -) d A n[g(A)}}, 



(21) 



The function q(lo) is defined via q(oj) = Q + {u) + g (lo) 



with Q^(cl>) = ± / dui' 



To keep the 



integrals convergent we introduce a high-energy cutoff, 
which in our case is provided by the charging-energy scale 
Eq. Then in the limit of Eq 3> max{w, k^T} we obtain 
for the real part of the £>-function, 



3? [g{oj)} = -2 aj { In 



fPE c 



V 2tt 



j3u 



(22) 



where ip denotes the digamma function. We identify the 
first term in Eq. (f2"Tj) as the contribution due to the renor- 
malization of the charging-energy gap A, 



(i,rcn) 



a 5ft [q{A)\ — (n 



0,o) 



(23) 



where the charging-energy gap is renormalized by A — > 
A + aoR[g(A)]. 

The broadening of the spectral function of the island 
due to the tunnel coupling to the leads yields a con- 
tribution to the first-order correction of the island oc- 
cupation number, which can be written as the sum of 
the contributions of the individual barriers, (71) ro = 

^(i.broad.L) + (n) (i,broad,R)_ The second term of Eq . QJJ 

can be identified as the contribution due to broadening, 



(i, broad) 



,0,o) 



i)ao0A»[e(A)] , (24) 



where d A denotes the derivative with respect to the 
charging-energy gap A. The amplitude of the broadening 
depends on the occupation of the island: if the island is 
unoccupied, the prefactor is — ^; if the island is occupied 
by a single electron, the coefficient is \. 



fa 1) 

The expressions of the matrix elements of kernels WJ ' ; 
and wf' , which are needed to compute the probabili- 
ties Eq. ([25]) . can be found in Appendix lAl Finally, the 
pumped current in first order in ao can be written as 



r(a,l) 



(i, broad, L) «0 



i (i,rcn) 



ao dt 



(26) 

The first term is related to the time variation of the cor- 
rection of the average occupation number due to broad- 
ening of the resonance at A induced by tunnel coupling 
the island to the left lead. It is, therefore, associated with 
tunneling processes where electrons tunnel through the 
left barrier. Equation (|26[) shows that this first contribu- 
tion to the current consists of the total time derivative of 
the correction of the average occupation number due to 
tunnel coupling to the left lead and thus the integral over 
a whole pumping period yields zero. The second contri- 
bution to the pumped current in first order in ao is due to 
the renormalization of the charging-energy gap and has 
the same form as the zeroth-order contribution to the 
current, Eq. (|20[) . with renormalized parameters. As the 
renormalized average occupation number of the island 
^0, ren ) depends on the charging-energy gap as well as 
on the tunnel-coupling strengths a^, we obtain a nonvan- 
ishing pumped charge for any choice of pumping param- 
eters. Particularly intriguing is the case of pumping with 

the two barriers, since in this case the current I^'°\t) 
vanishes and the first nonvanishing contribution to the 
current is due to the renormalization of the charging- 
energy gap, and therefore it is an effect purely due to 
quantum fluctuations and interaction. 



C. Pumped charge in the weak-pumping limit 

Now we calculate the pumped charge per cycle in the 
weak-pumping limit. We write the pumping parameters 
as Xi(t) = Xi + SXi(t), where Xi is the time average 
and SXi(t) the oscillating part. The weak-pumping limit 
corresponds to computing the current in bilinear order in 
the time- varying parts of the parameters. 
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1. Pumping with the charging- energy gap and one 
tunnel- coupling strength 

We start by considering one tunnel- coupling strength 
«o an d the charging-encrgy gap A as pumping parame- 
ters. Expanding the current 1^ (t) up to bilinear order in 
the time variation of the pumping parameters and inte- 
grating over one cycle, we obtain for the pumped charge 



Q 



(a) 



dpA 



\(i,0) 



a p [g(A)] 



d(3A 
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zeroth order 
PE C = 10 

PE C = 50 

(3 E c = 100 



(27) 
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where 771 = <5cvq ^Jp dt is the area of the pump- 



ing cycle in parameter space, 



.0) 



the instantaneous 



average occupation number, where the charging-encrgy 
gap has been replaced with its time average A, and 

«o = + "0 ■ 

The first contribution to the pumped charge arises 

from the current in zeroth order in ao, the second from 
the first-order correction of the current due to the renor- 
malization of the charging-energy gap as discussed in 
Sec. IIVBI Hence, the first term is the dominant one. 

In Fig. [3] the pumped charge is plotted as a function 
of the time average of the charging-encrgy gap for dif- 
ferent values of the high-energy cutoff. The pumped 
charge is symmetric around zero and has its maximum 
there. This peak is associated with transitions between 
charge states where the island is unoccupied or occu- 
pied by a single electron. The pumped charge vanishes 
for charging energies much bigger than the temperature 
(|A| ^> k-sT) without changing its sign. The overall 
sign of the pumped charge is determined by the traver- 
sal direction of the pumping cycle in parameter space. 
In Fig. [3] we can clearly see that the maximum of the 
pumped charge decreases for increasing high-energy cut- 
off due to the correction introduced by the renormaliza- 
tion of the charging-encrgy gap. 



2. Pumping with the tunnel- coupling strengths 

When we choose both tunnel-coupling strengths ag 
and «q as pumping parameters the contribution in zeroth 
order in cvrj vanishes, as already mentioned. We therefore 
obtain for the pumped charge in the weak-pumping limit 
only the contribution due to the renormalization of the 
charging-encrgy gap. 



V2 



v(i,0) 



(28) 



where 772 



c'T 



5a* 



d (5an 



dt is again the area of the 



Jo _ ""0 dt 

pumping cycle in parameter space. The pumped charge 
as a function of the time-averaged charging-encrgy gap is 



FIG. 3: (Color online) Pumped charge through the single- 
island system up to first order in ao in units of ePr/ictQ / Oq as 
a function of the time average of the charging-energy gap for 
different values of Ec- The pumping parameters are ao and 
A. The tunnel-coupling strength is oto — 0.01. 
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FIG. 4: (Color online) Pumped charge through the single- 
island system up to first order in ao in units of er^/ao as a 
function of the time average of the charging-energy gap for 
different values of Ec- The pumping parameters are q?o and 



plotted in Fig. [U It is zero for A = and antisymmetric 
around zero. For charging energies much larger than the 
temperature (|A| ^> ksT) the pumped charge vanishes. 
Increasing the high-energy cutoff Ec, the amplitude of 
the pumped charge is increased. Finally, we emphasize 
that by adiabatic pumping with the two tunnel-coupling 
strengths we can explore the effects of renormalization of 
the charging-energy gap A. The pumping mechanism is 
peristaltic but it occurs via renormalization effects. The 
sign change of the pumped charge around the resonance 
clearly distinguishes this mechanism from lowest-order 
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pumping with one barrier and gate voltage. We remark 
that such a sign change has been experimentally observed 
in charge pumping through carbon nanotubes with the 
help of surface acoustic wavesi^i There, this feature has 
been understood as adiabatic pumping through an effec- 
tive few-level double quantum dot with changing level 
positions. 
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ADIABATIC PUMPING THROUGH A 
DOUBLE-ISLAND SYSTEM 



FIG. 5: Sketch of the double-island system. It consists of 
two metallic islands tunnel coupled in series to noninteracting 
leads which are kept at the same chemical potential. We 
consider as pumping parameters the two charging-energy gaps 



In this section we consider pumping with the two 
charging-energy gaps in a system consisting of two 
tunnel-coupled metallic islands with no applied bias volt- 
age. This system corresponds to the single-electron pump 
of Pothier et al~ The left and right islands are labeled 
by L and R, respectively. Notice that we will use the 
same labels for the the left and right leads. Whenever 
confusion could arise, we specify whether we refer to an 
island or a lead. The charging energy of the system reads 



E ch (N h ,N R ,N x , L (t),N x , R {t)) 
-- E C l(N l -N x , l ) 2 + E C r(N r - 
+ E cm (Nl - N Xth )(N R - N X , R ) 



N 



x,R) 



(29) 



where Nurs is the total charge on the left (right) island 
and A^ l(r) the gate-induced charge on the left (right) 
island. The prefactors Eql, Eq R , and Ecm can be eas- 
ily computed as a function of the junction and gate ca- 
pacitances. At low temperature, the Hilbcrt space for 
the double-island system can be truncated to the three 
states |0, 0), |1, 0), and |0, 1), corresponding to no excess 
charge in the system, one excess charge on the left island 
and none on the right one, and one excess charge on the 
right island and none on the left one, respectively. The 
corresponding energies are Eofl(t), Eifiif), and Eq^i{€). 
We define the two time-dependent charging-energy gaps 
A L (t) = E lfi (t)-E ,o(t) and A R (t) = E ,i(t)-E 0fi (t) of 
the islands and the abbreviation Am(£) = Ar,(£) — Al(£). 
The occupation probabilities for these charge states are 
Po, Pl, and p R and are collected in the probability vector 
P = (po,Pl,Pr) T ■ A sketch of the double-island system 
is shown in Fig. [5] 

In comparison to the single-island system, the tunnel- 
ing events in the double-island system include additional 
tunneling between the two islands. We define the direc- 
tion of these tunneling lines to correspond to an electron 
tunneling from the right to the left island. Since no trans- 
port voltage is applied the pumped current vanishes in 
instantaneous order and we need to compute the first 
adiabatic correction. As in this system the current in 
zcroth order in the tunnel-coupling strength is the dom- 
inant term, wc restrict ourselves to this contribution. 



A. Pumped current in zeroth-order in ao 



To compute the pumped current in zeroth order in ao 

(a-l) 



we need the probabilities pj ' , determined by Eq. |T 
together with Eq. ||5J). For this purpose, only the instan- 
taneous kernels in first order in ao are necessary. Ex- 
amples of diagrams contributing to the matrix elements 
of the kernel Wj 1 ' 1 ' for the double-island system can be 
found in Appendix [B] Following this procedure we find 
for the instantaneous occupation probabilities 



(ho) 
Pt 



i 



= -/3A L + e -/3A R + I 



( 



\ 



1 

,-/8A R 



(30) 



Making use of the notation a* = a*(A s ) and a 
a,T with r € {L, R, M} as well as J\f = a^a 



a L a R -t-a L a R 



-a M a R 



' a M a R 



M 



= a+- 



for the denominator, the 



adiabatic probabilities can be written as 



(a,-l) 



2nAf 



«m + & R «M + a i 

a R + a~l 
Ql + «m a L + a R 



lit* (; o) 

a R + a M a R + aj | — p t ' 



(31) 



By substituting these results in Eq. (T3"|) we obtain the 
adiabatically pumped current in the left lead in zeroth 
order in ao ; 

where the functions and g R are defined by 



9b 



a L a M + Q?l 

N 

Ql«m + ajaj, + aLa R + a^a R 



(33a) 
(33b) 



The current, Eq. (|32[) . depends on the time derivative of 
the occupation probabilities of both islands. The current 
in the right lead can be obtained by exchanging L and R 
in Eq. (|32[) ; in doing so we have to replace also aj^ with 

a M- 
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B. Pumped charge 



To compute the pumped charge per cycle we simply 
need to integrate the current I^"°'(t), given in Eq. (f3"2| . 
over one period. Using simple mathematical manipula- 
tions, the pumped charge per period can be cast in the 
form 



Q 



Qal.Ar = e / dt 



dD L dAj 
dA~ dt 



where the emissivity for the left lead is given by 



dA~ 



9B 



dAi 



(i,0) 



9h 



dA, 



Following Ref. [l| the pumped charge per period can be 
written as an integral over the area A spanned by the 
cycle in parameter space and reads 



f ,A ,A f 9 dD 

Q^An=eJ A dA R dA h i [M -^ K 



d dD L 
9A L dA R 



(36) 

We now consider the weak-pumping regime and we write 
the pumping parameters, here the charging-encrgy gaps, 
as their mean values plus a small time-dependent varia- 
tion around it, i.e., A r (t) = A r + 5A r (t). 

Equation (|36[) in the weak-pumping limit yields for the 
zcroth-ordcr contribution to the pumped charge 



Qa l ,Ar 



j={R,L} 







Hi 







d Ar dA L 

dgj d 



d Al dA R 



(37) 



where % = J T SA R d s £ h dt is the area enclosed by the 
pumping cycle in parameter space. In the functions g~j 

and the instantaneous occupation probabilities (pj 1 ' ^- 
the charging-energy gaps have been replaced with their 
time averages Aj. 

A density plot of the pumped charge as a function of 
the time averages of the pumping parameters Al and 
Ar is shown in Fig. [S] It consists in a peak with its 
maximum in the origin (for symmetrically chosen tun- 



nel couplings Q!q 



«Q ). This can be seen more 



clearly in Fig. [7] where we show a cross section of the peak 
along the line with Al = — Ar (solid line) and along 
the line with Al = Ar (dashed line). Notice that the 
peak is symmetric for the case of antisymmetrically cho- 
sen charging-energy gaps, while the symmetry is absent 
for symmetrically chosen charging-energy gaps. Going 
around the origin of the plot in Fig. [6] with constant ra- 
dius R = \J A^ + Ap, the pumped charge exhibits a 
maximum when the polar angle takes the values j and 
|tt, which correspond to Al = Ar . 



A L , A R 



[e(3 2 ri 3 ] 



0.040 



(34) (3A R 



(35) 




FIG. 6: (Color online) Pumped charge through the double- 
island system in the weak-pumping limit in zeroth order in 
qo in units of e(3 2 rj3 as a function of the time averages of the 
pumping parameters Al and Ar in units of the inverse tem- 
perature = ksT. The white square marks the area in pa- 
rameter space over which we integrate to obtain the pumped 
charge in the strong-pumping limit. 



So far, we discussed the pumped charge in the weak- 
pumping limit. Now, we want to describe the transition 
from weak to strong pumping. Therefore, we increase the 
variation of the pumping parameters 5Aj(t) and, hence, 
the area enclosed in the parameter space. The pumped 
charge over a large pumping cycle can be obtained by 
integrating the weak-pumped charge of Eq. (|37[) over the 
area of the cycle. For the sake of definitcness, we choose 
as a pumping cycle a square loop around the origin, as 
shown in Fig. [6l The result for the charge as a function of 
the side of the loop, \f~A, obtained in this way is plotted 
in Fig. [5] In this figure we see that we obtain charge e if 
the size of the cycle is larger than k R T . This agrees with 
the results of the experiment of Pothier et al^ where 
charge quantization was measured for strong adiabatic 
pumping. Finally, we wish to stress that by means of 
our technique we are able to describe the transition from 
weak to strong pumping, recovering charge quantization 
for strong pumping cycles encircling the degeneracy point 
A L = A R = 0. 



VI. CONCLUSIONS 

We developed a diagrammatic real-time approach to 
adiabatic pumping through a system of tunnel-coupled 
metallic islands, performing a systematic perturbative 
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FIG. 7: Pumped charge through the double-island system in 
the weak-pumping limit in zeroth order in ao in units of ef3 2 ri3 
for Al equal to Ar (dashed line) and Al equal to — Ar (solid 
line) . 




FIG. 8: Pumped charge per cycle as a function of the side of 
the square cycle in parameter space obtained by means of a 
numerical method. 



expansion in powers of the tunnel- coupling strengths. 
This method allowed us to identify the different phys- 
ical processes which contribute to the pumped charge. 
We first applied our formalism to a single-island system 
consisting of one metallic island with Coulomb interac- 
tion tunnel coupled to two nonintcracting leads. We 
computed the pumped charge up to first order in the 
tunnel-coupling strength, finding that the contribution in 
first order is due to the renormalization of the charging- 
energy gap. For the case of pumping with the two tunnel- 
coupling strengths, this term becomes the dominant one. 
We emphasize that the single-island system works like a 
peristaltic pump if one tunnel-coupling strength and the 



charging-encrgy gap are the pumping parameters. For 
pumping with both tunnel-coupling strengths the peri- 
stalting mechanism occurs by means of charging-cnergy 
gap renormalization, therefore enabling experimental ac- 
cess to renormalization effects. 

Furthermore, we considered pumping with the 
charging-energy gaps in a system consisting of two 
tunnel-coupled metallic islands. We calculated the 
pumped charge through this double-island system and 
studied the issue of pumped-charge quantization. First, 
we investigated the pumped charge in the weak-pumping 
limit. For example, we found that in experiments a sym- 
metric choice of the average of the charging-energy gaps 
should be preferred to maximize the pumped charge. 
By integrating the weak-pumped charge we obtained the 
pumped charge over a larger pumping cycle. Therefore, 
we were able to describe the transition from weak to 
strong pumping. In the strong-pumping limit we found 
the charge to be quantized. This is consistent with the 
results obtained by Pothier et al^- for strong pumping 
with the single-electron pump. 

We acknowledge financial support from DFG via SPP 
1285. 



APPENDIX A: EXAMPLES OF DIAGRAMS FOR 
THE SINGLE-ISLAND SYSTEM 

In this appendix, we give the results for the matrix 
elements (Wt)i,j for the single- island system. These ma- 
trices are needed to compute the occupation probabil- 
ities of the metallic island. We start with calculating 
the matrix elements (Wj 1 )ij of the instantaneous ker- 
nels in first order in ao. The corresponding diagrams are 
shown in Fig. [91 where the additional dotted frequency 
lines must not be taken into account. As an example, 
we compute the two diagrams contributing to the matrix 
element (w{'' )o,i, shown in Fig. [Hb). The sum of the 



two diagrams reads 



(wP) 

V / o,i 



-23 



a r {uj) 



A - uj + i0 + 
= 2na-(A) . (Al) 

The instantaneous kernel in first order in ao reads 



W 



(U) 



2tt 



-o+(A) a" (A) 
a+(A) -a" (A) 



For the derivative of the instantaneous kernel in first or- 
der in a o we get 



3WJ U) =a d A n[g(A)} 



1 -1 
-1 1 
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FIG. 9: (a) (wf /a,1) )o,o , (b) (Wf /a,1) ) ,i , (c) (Wf /a ' 1} )i,o , and (d) (W t (i/a,1) )i,i . Diagrams contributing to the instanta- 
neous/ adiabatic matrix elements in first order in ao. The dotted frequency lines are only needed for the adiabatic corrections 



To compute the adiabatic correction we need to introduce 
additional external frequency lines according to rules (7a) 
and (8a) , which are drawn as dotted lines in Fig. [5J The 
evaluation of these diagrams leads us to the adiabatic 
correction to kernel in first order in ao, 



/ 



V 



1 -1 



Finally, we calculate the second-order contribution in 
cv to the matrix elements of the instantaneous kernel. 
As an example, we show all diagrams contributing to the 
matrix element (wf ,2 ')o,i in Fig. [lOl We have to sum 
over all indices r, r' and to integrate over all frequencies 
u>,u>'. Performing the limit z = 0+ and using Cauchy's 
principal value, we get 



(wf< 2 )) 



0,1 



-Trao a(A) d A 5i[£»(A)] 
+27ra »[g(A)] d A a" (A) 



For the matrix elements of the instantaneous kernel in 
second order in ao we find 



W 



(i,2) 



7ra a(A) d A U[g(A)} 



{ 



\ 




+27ra »[e(A)] <9 A 



-a+(A) a - (A) 
a+(A) -o-(A) 



APPENDIX B: EXAMPLES OF DIAGRAMS FOR 
THE DOUBLE-ISLAND SYSTEM 



In this appendix we give the results for the matrix el- 
ements (Wt')i j of the instantaneous kernel in first or- 
der in ao for the double-island system. Examples of the 
corresponding diagrams are shown in Fig. 1111 Their eval- 
uation is similar to the one for the instantaneous matrix 
elements in first order in ao for the single-island system, 
see Eq. ([All . 



V / 10,01 



= -23? 



dto- 



«m( w ) 



uj + iQj 



2t%(Am) 



Proceeding in the same way for the other matrix el- 
ements and using the notation af = a^(A s ) with s = 
{L, R, M}, we obtain for the matrix for the instantaneous 
kernel in first order in ao, 



W 



0,1) 



2tt 



— at — at 



a 



- a M - a L 



M 
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a M 
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